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It has been pointed out that the null energy condition can be violated stably in some non-canonical 
scalar-field theories. This allows us to consider the Galilean Genesis scenario in which the universe 
starts expanding from Minkowski spacetime and hence is free from the initial singularity. We use 
this scenario to study the early-time completion of inflation, pushing forward the recent idea of 
Pirtskhalava et al. We present a generic form of the Lagrangian governing the background and 
perturbation dynamics in the Genesis phase, the subsequent inflationary phase, and the graceful 
exit from inflation, as opposed to employing the effective field theory approach. Our Lagrangian 
belongs to a more general class of scalar-tensor theories than the Horndeski theory and Gleyzes- 
Langlois-Piazza-Vernizzi generalization, but still has the same number of the propagating degrees of 
freedom, and thus can avoid Ostrogradski instabilities. We investigate the generation and evolution 
of primordial perturbations in this scenario and show that one can indeed construct a stable model 
of inflation preceded by (generalized) Galilean Genesis. 

PACS numbers: 98.80.Cq 


I. INTRODUCTION 

Inflation in the early Universe [1, 2] is now an indis¬ 
pensable ingredient of modern cosmology not only to ex¬ 
plain the global properties of homogeneous and isotropic 
space with a vanishingly small spatial curvature but 
also to account for the origin of the primordial curva¬ 
ture perturbation that seeded cosmic structure forma¬ 
tion [3]. At present, despite the significant progress in 
the state-of-the-art precise measurements of the cosmic 
microwave background radiation (CMB) by WMAP [4, 5] 
and Planck [6, 7] missions, there is no single observational 
result in conflict with the single-field inflation paradigm 
[2]. In particular, the anti-correlation of the temperature 
and the E-mode polarization anisotropies on large scales 
observed by the WMAP mission strongly supports the 
superhorizon perturbations suggested by inflation [8]. 

In other words, once inflation sets in, virtually all the 
available cosmological observation data can be explained 
simultaneously irrespective of the initial condition of the 
Universe. This does not mean that we may be indif¬ 
ferent to the initial condition of the Universe before in¬ 
flation. On the contrary, in order to achieve complete 
understanding of the cosmic history, we must work out 
the very beginning of the Universe that may smoothly 
evolve into the inflationary phase. 

As is well known, as long as the null energy condi¬ 
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tion (NEC) is satisfied in the expanding phase, the Hub¬ 
ble parameter and the energy density of the universe in¬ 
crease backward in cosmic time. So, it is often claimed 
that, if one tries to discuss what happened before infla¬ 
tion and/or how inflation started, one needs to know the 
information of very high energy physics, and challenge 
the initial singularity problem [9] in terms of quantum 
gravity. But, this is not always the case. 

Recently, it was recognized that, if an action includes 
higher derivative terms of a scalar field like the Galileon 
terms, the NEC can be violated without encountering 
ghost nor gradient instabilities. See, e.g.. Ref. [10] for a 
recent review and Ref. [11] for a subtle issue of nonlinear 
instabilities. If the NEC is violated, the energy density 
can grow as time proceeds, contrary to the conventional 
wisdom. In the NEC violating theories, the universe can 
therefore start from the static zero-energy state described 
by the Minkowski spacetime from infinite past [12], and 
the universe starts expansion with the increase of the 
energy density. 

Such a picture of the emergence of the universe was 
first proposed by Creminelli et al. [13] with the name 
Galilean Genesis. In their model, however, the hot big 
bang state was postulated to be realized after the effec¬ 
tive field theory description breaks down as the energy 
density blows up beyond its realm of validity. Therefore, 
the theory to describe the most important epoch of the 
early universe is lacking there. 

Nevertheless, since their original idea is so interesting 
that a number of extension has been made in a wider 
class of scalar field theories [14-18] and various aspects of 
the Genesis scenario have been explored in the literature 
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[19-23], such as avoidance of the superluminal propaga¬ 
tion of perturbations and absence of primordial tensor 
perturbations. They have been unsuccessful, however, to 
realize transition from the Genesis phase to the hot big 
bang state within their model Lagrangians. 

In this paper, we take a different approach, namely, to 
make use of the Galilean Genesis to explain the initial 
condition of the Universe before inflation and smoothly 
connect it to the inflationary phase, thereby solving the 
initial singularity problem [9] and the trans-Planckian 
problem [24] (see also [25]) in inflationary cosmology. 

In fact, such an approach has also been put forward 
by Pirtskhalava et al. [26] recently. Their model La- 
grangian, however, gives rise to gradient instability as 
it is, although it has been argued there that higher-order 
structure of the effective field theory for perturbations 
possesses enough freedom to cure the gradient instabil¬ 
ity. Discussion on termination of inflation and reheating 
is not presernted there, either. 

In the present paper, we construct a specific model free 
from any catastrophic instabilities and with subluminal 
velocities of primordial perturbations. In our setup the 
universe starts from the Minkowski spacetime from infi¬ 
nite past and is smoothly connected to the inflationary 
phase followed by the graceful exit. For this purpose, we 
provide a generic Lagrangian capable of describing the 
background and perturbation evolution in all the above 
phases instead of choosing the effective field theory ap¬ 
proach because the latter cannot capture the evolution of 
the background and perturbations from pre-inflationary 
Genesis to the exit from inflation with the same single 
Lagrangian. 

Although we start with asymptotically Minkowski 
space at the past infinity for aesthetic beauty, it has been 
shown that the Galilean Genesis solution is an attractor 
for a variety of initial conditions including those with a 
negative Hubble parameter and/or finite curvature, pro¬ 
vided that the time derivative of the scalar field has the 
right sign [18]. 

The Horndeski theory [27] or the generalized Galileon 
[28], whose mutual equivalence was first shown in [29], is 
known to be the most general scalar-tensor tensor theory 
with the second-order field equations, and thereby avoid 
Ostrogradski instabilities in spite of having higher deriva¬ 
tive terms in the action. The theory can be generalized 
to have second-order field equations only in a specific 
gauge while maintaining the number of propagating de¬ 
grees of freedom. This possibility was realized recently 
by Gleyzes et al. [30] (see also Ref. [31]) and was ex¬ 
tended further by Gao [32]. The number of propagating 
degrees of freedom in these theories is indeed shown to be 
the same as that of the Horndeski theory [30, 32-36]. In 
this paper, we use the subclass of Gao’s framework as a 
concrete realization of the unified scenario starting from 
Galilean Genesis through inflation to the graceful exit. 

This paper is organized as follows. In the next section, 
we give a framework of our model and derive the back¬ 
ground equations of motion and the quadratic actions 


of cosmological perturbations. In Sec. HI, a concrete 
Lagrangian is constructed to describe our scenario be¬ 
ginning from the Genesis phase through the inflationary 
one to the graceful exit, and such a background dynamics 
is presented explicitly. In Sec. IV, we discuss the stabil¬ 
ity during each phase based on the quadratic actions of 
cosmological perturbations. In Sec. V, a concrete realiza¬ 
tion of our scenario is given. The final section is devoted 
to our conclusions and discussion. 


II. GENERAL FRAMEWORK 

Let us start with describing the general framework to 
construct and study our explicit realization of the early- 
time completion of inflation. We would like to consider 
theories composed of a metric and a single scalar field 
/i, and hence it will be appropriate to work in the Horn¬ 
deski theory. The Lagrangian of the Horndeski theory is 
of the form 

C = V^[G2{^,X) - G3{^,X)a^ + G4{^,X)R^^'> 
+G 5 (</,A)gWv'^V> + ---], (I) 

where X := is the four-dimensional 

Ricci scalar, and G^D is the four-dimensional Einstein 
tensor. We have four arbitrary functions of (j) and X 
in the Horndeski theory. This is the most general La¬ 
grangian having second-order held equations. Neverthe¬ 
less, it will turn out that this framework is insufficient 
for our purpose, and hence we have to go beyond the 
Horndeski theory. 

One can generalize the Horndeski theory to possess 
higher order held equations while maintaining the num¬ 
ber of propagating degrees of freedom [30]. The hrst step 
to do so is to perform an ADM decomposition by taking 
(j) = const hypersurfaces as constant time hypersurfaces. 
In the ADM language, the metric is written as 

ds^ = + 7y (dx* -b N^dt) {dx^ + N^dt) . (2) 

By dehnition is a function of only t, (j) = 4)(t), and 
X = (jp'/2N‘^^ where a dot denotes differentiation with 
respect to t, so any function of (j) and X can be regarded 
as a function of t and the lapse function N, provided 
that <j> and N~^ never vanish. Then, the Horndeski La¬ 
grangian (I) can be written in terms of the ADM vari¬ 
ables as £ = '^g^La with 

L2=A2{t,N), L3 = A3{t,N)K, 

U = A4(t, N) (K^ - Kl) + B4{t, N)R, 

£5 = A5(t, N) {K^ - iKKl + 2K%) 

+H5(£ N)K^^ (^Rij - , (3) 

where K^j and Rij are the extrinsic and intrinsic curva¬ 
ture tensors on the constant time hypersurfaces, and A 4 , 
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A^, B 4 , and are subject to the relations 


A 4 


-B 4 - N 


dB4 


A 5 


NdB 5 


( 4 ) 


Variation of the above Lagrangian with respect to N 
gives a second-class constraint that eliminates only one 
degree of freedom, as opposed to general relativity. The 
key trick to generalize the Horndeski theory is to notice 
that this property remains the same even if one liberates 
A 4 and A 5 from the restriction imposed by Eq. (4) [30]. 
We thus arrive at the so called GLPV theory that is 
more general than Horndeski but has the same num¬ 
ber of propagating degrees of freedom. One can move 
back to a covariant form of the Lagrangian by introduc¬ 
ing the unit normal to the constant time hypersurfaces as 
= —d^(j)l^/ 2 X^ writing the extrinsic curvature tensor 
in terms of n^, and using the Gauss-Codazzi equations. 
Since there are six arbitrary functions of t and N in the 
ADM form, the resultant covariant Lagrangian has six 
arbitrary functions of </> and X. 

The above idea has been pushed forward by Gao [32], 
who proposed a unified framework to study single scalar- 
tensor theories beyond Horndeski. One can write a gen¬ 
eral Lagrangian in the ADM form as 


B — Y^iV [c?o T diR -f d 2 R^ -I- • ■ • -I- (ct-o T ctiR -I- ■■■) X 
+ (a2X^ +---)K^j+ biK^ + b2K^jK^^ +•••], (5) 


where the coefficients doi di, ... are arbitrary functions 
of t and N. The Hamiltonian depends nonlinearly on 
N as in the GLPV theory, giving rise to a single scalar 
degree of freedom on top of the traceless and transverse 
gravitons [35]. 

In this paper, we will employ the Lagrangian C = 
with 

L2 = A2it,N), L3=A3{t,N)K, 

L 4 = A4it,N) - Kl) + B4{t,N)R, 

L 5 = A,{t, N) {\2K^ - iX^iKKl + 2Kl) 

+B^{t, N)K^^ [r,j - 1 ( 7 ,,, ( 6 ) 


where Ai, A 2 , and A 3 are constant parameters of the the¬ 
ory. This is a deformation of the GLPV Lagrangian and 
belongs to a subclass of Gao’s framework. The gener¬ 
alization to this level is sufficient for the purpose of the 
present paper. The GLPV theory is recovered by taking 
Ai = A 2 = A 3 = 1. 

Given the Lagrangian ( 6 ) in the ADM form, one can 
restore the scalar degree of freedom (j) to write its covari¬ 
ant expression in the same way as in the GLPV theory. 
However, it will be more convenient for our purpose to 
use the explicitly time-dependent Lagrangian, because 
by doing so one can easily design the Lagrangian so as 
to admit the desired cosmological evolution. 

Before specifying the suitable form of A 2 {t,N), 
A^{t, N), ... to construct our early universe model, let us 


derive the general equations governing the background 
and perturbation dynamics of cosmologies based on the 
Lagrangian ( 6 ). The ADM variables are given by 

N = N{t) (1 + Sn ), Ni = NdiX, 

Xij — d ^ j (7) 

where C is the curvature perturbation in the nnitary 
gauge and hij is the transverse and traceless tensor per¬ 
turbation. A spatially flat background has been assumed 
and the spatial diffeomorphism invariance was used to 
write 7 y in the above form. In the following, the back¬ 
ground value of the lapse function is denoted by N where 
there is no worry about confusion. 


A. Background Equations 

Substituting Eq. (7) to the Lagrangian ( 6 ), we obtain 
the background part of the Lagrangian as 

£(o) = (A 2 -f 3 A 3 H-f 6774 ^ 477 ^-f 6775 ^ 577 ^) ( 8 ) 

where 774 := (3Ai — l)/2, 775 := ( 9 A 2 — 9 A 3 -I- 2)/2, and 
77 := d/{Na). At the background level, Ai, A 2 , and 
A 3 just rescale A 4 and A 5 . In what follows we simply 
consider the case with 774 > 0 Ai > 1/3. Since we 
are considering a spatially flat universe, we have Rij = 0 
at zeroth order, and hence B 4 and B^ play no role in the 
background dynamics. Varying Eq. ( 8 ) with respect to 
N and a, we obtain, respectively, 

-£ ■.= {NA 2 )'+ ?>NXH+ Qri 4 N^{N-^A 4 )'H‘^ 

+6r]5N^{N-^A5yH^ 

= 0, (9) 

r := A 2 - 6774^4772 - 12775 A 5773 

{As + 4774A477 -t erjsAsH'^) 

= 0 , ( 10 ) 

where a prime represents differentiation with respect to 

TV. The background equations contain at most second 

derivatives of the scale factor and first derivatives of the 
Lapse function. 


B. Cosmological Perturbations 


The quadratic Lagrangian for the tensor perturbation 
is given by 


jriA 


Na^ 


Gt 12 ,2 


where 


( 11 ) 


Gt 

:= - 2 A 4 - 6 ( 3 A 3 -2)^577, 

( 12 ) 


1 dSs 

(13) 

Tt 

- 2 ^^+TVdf 
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The equation of motion contains at most second deriva¬ 
tives both in time and space. The tensor perturbation is 
stable provided that Gt > 0 and Jr > 0 . 

The quadratic Lagrangian for the scalar perturbations 
is given by 


/’(2) 


Na'^ 


-3Ga 


a2 


{dCf + 


-2eSn^ + 2gAiT^+ GQSn ^ 


a? N a? 

-2gB5n^-C^^ 


N 


(14) 


where the coefficients are defined as 


s 

:= VAo-f-h 

22 -^2 



-f3r/4 ( 2 A 4 - 2NA'^ + N‘^A'1) 

+ 3 r ;5 ( 6 A 5 - 4VA^ + V"A") H^, 

(15) 

0 

N A' 

■■= ^ - 2774 (A 4 - VA;) iJ 



-3r75(2A5-iVA^)if", 

(16) 

Ga 

:= -2774^4 - erj^A^H, 

(17) 

Gb 

:= 2 {B4 + NB'^)-HNB'^, 

(18) 

C 

:= (1-Ai)A4-(6 + 9A2-15A3)A5H, 

(19) 


We require that t/s > 0 in order to avoid ghost insta¬ 
bilities. However, we allow for a negative sound speed 
squared, cl := Fs/Gs < 0 , for a short period of time. 
In the absence of the term (C = 0), a negative sound 
speed squared would cause a rapid growth of instabilities 
for large k modes. In this paper, we consider theories 
with C 0 , so that the curvature perturbation with large 
k can be stabilized by requiring that 'H.sjGs > 0 . 

As will be seen in the rest of the paper, the sound 
speed squared becomes negative at the transition from 
one phase to another. Such a behavior should not occur 
even for a tiny period because high wavenumber modes 
would grow exponentially rapidly. However, we could not 
avoid it not only within the Horndesky theory but also 
the GLPV theory despite we analyzed extensive models. 
On the other hand, we have not been successful in prov¬ 
ing that this is an inevitable consequence. Since our pri¬ 
mary purpose is to show an existence proof of the model 
to realize our intended cosmic evolution without any in¬ 
stabilities, we construct a specific model by going beyond 
the GLPV theory and invoking the k^ term. 


III. STARTING INFLATION FROM 
MINKOWSKI 


A. Construction of the Lagrangian 


and note the relation Gt = Ga — 3C. One has C = 0 in 
the Horndeski and GLPV theories, in which Ai = A 2 = 
A 3 = I. Therefore, the last term in the Lagrangian (14) 
is the novel consequence of theories beyond GLPV. 


From 6 C 

f^/S{Sn) = 

6 n = 

1 

02 -hSC 

d^X 

1 

a" 

02 -hSC 


C 

N 


( 30 "+ £ 05 ) 4 - 0 ^ 5 ^ 

N 


, ( 20 ) 
( 21 ) 


Substituting Eqs. (20) and (21) into Eq. (14), we obtain 
the reduced Lagrangian for the curvature perturbation. 


c(2) 


NF 


Ssj^ + C[:Fs^ 



C , 


( 22 ) 


where 


Gs 

,oc 

■" 02 + SC+^^^’ 

(23) 

Fs 

Id/ ciQGbGt^ 

(24) 

Nadt iv02-hEC ) 

Fs 

QIC 

■ 02 -hSC' 

(25) 


Thus, if C 0, the equation of motion for C has the 
fourth derivative in space, giving the dispersion relation 


w 


2 


J 2 , 

Gs Gsa^' 


(26) 


The Lagrangian we study in this paper is character¬ 
ized by a single time-dependent function f(t) and four 
functions 02 , 03 , 04,05 of N: 


A 2 


(27) 

A 3 

= M|/-("“+')a3(iV), 

(28) 

A 4 


(29) 

A 5 

= a^{N), 

(30) 


where a (> 0) is a constant parameter. We have intro¬ 
duced the mass scales Ma (and the Planck mass Mpi), 
so that f{t) and aa{N) are dimensionless. The other two 
functions, B 4 and B^, are arbitrary at this stage because 
they have no impacts on the background dynamics. Note 
that / is not a dynamical variable. Specifying the func¬ 
tions / = f{t) and Oa = Oa(N) amounts to defining a con¬ 
crete theory. The above forms of Aa are chosen so that 
the theory admits an inflationary universe preceded by 
the generalized Galilean Genesis while retaining much of 
the generality. Other choices could be possible and hence 
we do not claim that this is the most general description 
of such scenarios at all. Instead, as we mentioned above, 
we would provide the existence proof of desired models 
by demonstrating that a sufficiently wide class of healthy 
models can indeed be constructed. 

We design f{t) so as to implement the (generalized) 
Galilean Genesis followed by inflation and a graceful exit 
from the prolonged inflationary phase. Our choice is 

/ ~ lot (/o = const < 0 ) 


(31) 
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well before t = to, and 

f - fi= const (32) 

for t ^ to- As our time variable starts at t = —oo with 
asymptotically Minkowski spacetime configuration, t is 
large and negative in the beginning, so we find / ^ 1 
in Eq. (31). As will be seen shortly, the initial stage 
described by Eq. (31) corresponds to the generalized 
Galilean Genesis, while the subsequent stage described 
by Eq. (32) to inflation. After a sufficiently long period 
of the inflationary stage, we assume that 

/ ^ tl/(a+l) (33) 

foe f ^ fend, where Gnd is the time at the end of inflation. 
With this the universe exits from inflation. In what fol¬ 
lows we will investigate the background evolution of each 
stage. 


B. Genesis Phase 

Assuming that H ~ in the first stage where 

/ is given by Eq. (31), let us look for a consistent solution 
for large /. The background held equations read 

-S = + =0, (34) 

= 0. (35) 

It can be seen from Eq. (34) that the lapse function N is 
a constant, N = Nq, satisfying 

a 2 (iVo) + Aoa' 2 (Ao) =0. (36) 

Then, H is consistently determined from Eq. (35), which 
can be written as 


27?4Mgj dH 2ia+l)^ ^ Q 
No dt 


(37) 


where 


P — Af 2 02 (^ 0 ) -I- (2a -I- l)M|a 3 (A’o)-^ (38) 

i\o 

is a constant. This leads to the generalized Galilean Gen¬ 
esis solution [18]: 


H 


a 


_ P _ Afp „_(2a+l) 

2(2a + l)r74M2j/o|^ 

_ P M f-2a 

4a(2a -I- l)r]4Mp^ 


1 

(-t)2“+l 


(39) 


(40) 


It is required that p/p 4 < 0 to guarantee H > 0. We have 
thus arrived at the generalized Galilean Genesis solution 
starting from the Lagrangian written in the ADM form 


rather than in the covariant form. The original Galilean 
Genesis solution found in Ref. [13] corresponds to a = 1. 

In deriving the above solution, Mlf~'^°‘a 4 (c A 4 ) and 
Mofao (= A 5 ) are always subdominant due to the as¬ 
sumed scalings ~ f~^°‘ and ^ /. Therefore, any choices 
of 04 (iV) and 05 (iV) will not spoil the above Galilean 
Genesis solution. As will be seen in the next section, 
those two terms are also irrelevant to the stability condi¬ 
tions during the Genesis phase. 


C. Inflationary Phase 

The Galilean Genesis phase will end at t ^ to since the 
function / is constant for t > to- In the subsequent phase 
we obtain the de Sitter solution, N = Nini = const and 
H = Hint = const, satisfying 

-8 = (Mnf A 2 )' + 3M„fA^/7i„f + dr^4Nli{NrjA4)'Hl^ 

+6rjoNUNrjAoyHlt = 0 , ( 41 ) 

V = A 2 - = 0. ( 42 ) 

(Note that Aa is now a function of N only and is inde¬ 
pendent of t.) 

A t-independent Lagrangian in the ADM form can be 
recast in a covariant Lagrangian with the shift symmetry, 
<f) ^ (j) + c- This implies that the above exact de Sitter 
solution corresponds to kinetically driven G-inflation. If 
one invokes a weak time-dependence in /, one obtains 
quasi-de Sitter inflation instead. 


D. Graceful Exit 

After the prolonged phase of inflation, / is given by 
Eq. (33). We assume that t is sufficiently large, so that 
/ ^ 1. Then, we have a consistent solution with N = 
Ne = const and ~ l/t^ ^ y- 2 (a-i-i) ^ satisfying 

-8 = (fVeA2)'-f 3r74M|ii72 + c>(/-(3«+2)) 

= 0, (43) 

V = A 2 + 3r?4M2ji72 + + 0(/-(3“+2)) 

= 0. (44) 

Thus, one can implement a graceful exit from inflation. 

It follows from Eq. (43) that 

{N,a2y < 0. (45) 

It can be shown using Eqs. (43) and (44) that, during 
this third stage. 


oc 


m := 


3Nea2 

{Nea 2 y ■ 


(46) 


It is therefore necessary to impose m > 0 < 0. 
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In the standard potential-driven inflation models [2] 
inflation is followed by coherent field oscillation of the 
inflaton scalar field which decays to radiation to reheat 
the universe. In the present approach the scalar field (p 
is used to specify constant time hypersurfaces, so that 0 
may not vanish in order to preserve one-to-one correspon¬ 
dence between 0 and the cosmic time t. Hence one must 
switch from the ADM language we used to construct the 
action to the conventional “0 language” at this point in 
order to apply the standard reheating mechanism, which 
is all right but looks like sewing a fox’s skin to the lion’s. 

Here instead we consider another reheating mechanism 
which can take place without breaking the one-to-one 
correspondence between 0 and t, namely, the gravita¬ 
tional reheating due to the change of geometry or the 
cosmic expansion law [37-41]. 

During the transition from the de Sitter inflation to 
a decelerated power-law expansion, conformally non¬ 
invariant particles are produced with the initial energy 
density 


IV. PRIMORDIAL FLUCTUATIONS AND 
STABILITY 

Having obtained the background evolution of our sce¬ 
nario, let us investigate the nature of primordial pertur¬ 
bations and stability, using the result of the generic anal¬ 
ysis in Sec. HB. 


A. Genesis Phase 

During the Genesis phase, we have 

]\/f3 

0 ~ ^/-(2«+i)iVoa' + 

C~^(Ai-l). (52) 


Pr = (47) 

where ct is a factor determined by the effective number of 
conformally noninvariant fields and the change of the ge¬ 
ometry. For example, for m = 6 or 4, a single minimally 
coupled massless scalar field contributes to a by 


"'“si?'"(us)’ = 


respectively [41, 42]. Here At is the time required for 
the transition. In case it is nonminimally coupled with 
a coupling parameter a factor (1 — 6^)^ is multiplied 
there. 

In order for the radiation thus created to dominate 
the universe, the energy density of the scalar field must 
dissipate more rapidly, namely, 

m > 4 402 + A[ga 2 > 0, (50) 

then, the reheating temperature at the radiation domi¬ 
nation is given by 


Tr 


f 30 y/4 

WgJ V 3 ) 



2 l(m-i) 

77inf 5 
(51) 


where g* is the effective number of relativistic degrees of 


freedom and we have assumed the universe would evolve 


in the same way as in the Einstein gravity after infla¬ 
tion. If long-lived massive particles are copiously pro¬ 
duced at the gravitational particle production, the re¬ 
heating temperature may be significantly higher then the 
above value. Furthermore, the decay of quasi-flat direc¬ 
tion may produce a large amount of entropy to reheat 
the universe efficiently and create matter particles [43]. 


Obviously, the kinetic term of the tensor perturbations 
has the right sign, Gt > 0. For large /, we see EC 0^ 
(as long as C ^ 0), and hence 

G"^ Gl 

Gs-^ + 5 Gt, -Hs ^ (53) 

This implies that Gs — const, while Hs ^ (— 

The kinetic term of the curvature perturbation has the 
right sign if 

O \ _ 1 

es >0 T^> 0 - (54) 

Ai — i 

Thus, it is sufficient to impose 

Ai > 1. (55) 

(We are considering only the case with Ai > 1/3.) An¬ 
other stability condition, "Hg > 0 , is equivalent to requir¬ 
ing that 

(N^a'^y > 0. (56) 

Since Tt depends on B 4 and H 5 and these two func¬ 
tions are irrelevant to the background dynamics, the con¬ 
dition Ht > 0 can easily be satisfied without spoiling the 
Genesis background. Suppose for simplicity that 

B4 = ^, B, = 0, (57) 

where /3 (> 0 ) is a constant. Then, Tt = Gb = > 

0. For the scalar perturbations we have 

o.a.2 [M|a2 + (2a + l)M|(/o/iVo)(iVoa3)' l' 

^ [ (2a + l)(Ai - l)M|(Y2a')' 2_ 

= const. (58) 
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FIG. 1: Schematic diagram of the behavior of curvature per¬ 
turbation in {y, a/k) plane with y decreasing toward the right. 
In the region below (above) the red broken curve, is domi¬ 
nated by the term proportional to (k^)- Modes with k < k, 
experience the break down of the WKB approximation around 
the point crossing the blue solid curve beyond which ( is 
frozen, while modes with fc > fc* do not. 


This can also be made positive by an appropriate choice 
of 03 (A^). It should be noted that if 03 = 0 then we 
inevitably have J^s < 0 ; the L 3 term is crucial for the 
stable violation of the NEC. Note also that, if we take 
sufficiently small /3, the sound speed Cs can be smaller 
than unity, which applies also to the other two phases 
discussed below. 

Let us move to discuss the nature of the primordial 
fluctuations in the Genesis phase. Since Gt ^ ^ 

const, the tensor perturbations behave in the same way as 
in the Minkowski spacetime. Therefore, no large tensor 
modes are generated during the first stage of our scenario. 

The behavior of the curvature perturbation turns out 
to be more nontrivial, as sketched in Fig. 1. Recalling 
that Gs ^ const, Ts ^ const, and 'H.s ^ (—the 
equation of motion for () in the Fourier space is of the 
form 


dy2 


+ — 0, 


(59) 


where y := —NqI > 0 and 


-I- ki°'k‘^y 


2oi 7,4-,2a+2 


(60) 


with Cs and k^ being some constants. For sufficiently 
large y, we have ~ kl°‘k‘^y^°‘~^^. One may define 
the time at which this approximation breaks down as 

2/break := , and for y < l/break 

we have ~ c^k^- 

With some manipulation, it is found that 




d^w/di/^ 



(61) 


where /c* := Cs This implies that for the modes 

with k > k^ the WKB approximation is always good in 
the Genesis phase, 

Cfe oc ^ 

Vw 

giving Cfc oc for y < //break- Thus, the am¬ 

plitude of those modes at late times in the Genesis phase 
is given by 


exp 


ujdy' 


(62) 




k^ 

GsCs 


{k > /c*). 


(63) 


For the modes with k < k^:, the WKB approxima¬ 
tion breaks down at some time and then the curvature 
perturbation freezes. This “horizon crossing” occurs at 
y ^ yfreeze '■= k ~‘^^It Can be seen that 

y freeze > 2/break for k < fc*,^ which allowS US tO Study 
the freezing process by using the solution to Eq. (59) 
with Ri fcJ“A:^2/^“~''^- Tde exact solution in this case 
that matches the positive frequency WKB solution for 
y » 2/freeze IS given by 

Cfc oc y^/‘^Hl^\-2vk‘ik^y°‘+‘^), v := - ^ , (64) 

2[a -I- 2) 


where is the Hankel function of the first kind. The 
frozen amplitude can thus be evaluated by taking the 
limit y ^ y freeze In the solution (64), leading to 


kW 


^2 / ^ \ (3a-|-4)/(a-|-2) 

Gs V^*/ 


(65) 


For y < ybreak, uj is dominated by the Csk term where 
the solution (64) is no longer exact. The frozen ampli¬ 
tude (65), however, is still valid even in this regime since 


^ Note in passing that ^break = 2/freeze = ^ for the k = 

mode. The Genesis phase could end sufficiently early so that 
—A^O^O > If this is the case, we only need to care 

about the modes with /c < /c*. 

















the solution to Eq. (59) with the effective frequency (60) 
does not oscillate any more and remains constant. Hence, 
the expression of the power spectrum (65) is correct for 
the entire range of A: < fc*. 

To summarize, the power spectrum of the curvature 
perturbation generated during the Genesis phase is blue 
and hence is suppressed on large scales. 

B. Inflationary Phase 


(4/7r)|r(5/4)|2a;-3/2, so that one can take the limit 
Cg —>• 0 smoothly to get 

ttH? 

^ 8es|r(5/4)|2e3/2- (72) 

We have approximated the inflationary phase as exact 
de Sitter. If we consider a background slightly different 
from de Sitter by incorporating weak time dependence in 
/, we would be able to obtain a tilted spectrum of C- 


In the (de Sitter) inflationary phase, Qt, Gs, 
and Hs are time-independent. We require that all those 
coefficients are positive during inflation in order to avoid 
instabilities. 

Since the quadratic action for the tensor perturbations 
is essentially the same as that of generalized G-inflation, 
the power spectrum of the primordial tensor perturba¬ 
tions is given by [29] 


T’t 


Gt^^ H? 


inf 


47r2 


( 66 ) 


The equation of motion for the canonically normalized 
variable Uk ■= V‘^Gs(iCk during inflation is of the form 


dr^ 




(67) 


where 


= c2fc2 -I- 


2l,4_2 


( 68 ) 


with Cg = Fs/Gs and e := IGy^ being dimen¬ 

sionless constants. Here, we have introduced the confor¬ 
mal time T (< 0) defined by adr = NAt. The dispersion 
relations of this form have been studied in the context of 
inflation, e.g., in Refs. [44, 45]. The positive frequency 
modes are given by [45] 


Uk 


(-2efc2T)i/2 


(69) 


C. Graceful Exit 


After inflation, we have Gt — Ft = 


Fs 

5'lf2 1 + W2 

^ Xi-i + e ’ 

(73) 


, ^0 3Ai - 1 

(74) 

Gs 


ns 

yMi Ai - 1 e 

m 3Ai - 1 Ai - 1 -f 

(75) 


where to simplify the expression we introduced 

4 (iVea 2 )' 




3{Ny^y 


(76) 


Recalling that we have been imposing Ai > I, all of these 
coefficients are positive provided that £m > 2(Ai — I). 
This condition can be written equivalently as 


Nea '2 



I) (< 0). 


(77) 


V. A CONCRETE EXAMPLE 


Let us provide a concrete Lagrangian exhibiting the 
Genesis-de Sitter transition. The Lagrangian is charac¬ 
terized by 


^ ^0 _ 7 

3A4 ’ fV3 ’ 


(78) 


where IT„_m is the Whittaker function. Taking the limit 
T —>■ 0 , the power spectrum of the curvature perturbation 
can be calculated as 




Hit 


2GsclF{cl/ey 


(70) 


where Nq (> 0) and 7 (> 0) are constants. We take 04 = 
05 = 0, H 4 = Ml/2, and H 5 = 0. We also take Ai > 1 to 
guarantee the stability. This corresponds to the (Ai > 1 
generalization of the) unitary gauge description of the 
Lagrangian considered in Ref. [13]. In the Genesis stage 
we have 


where 


F(x) := lx- 3 / 2 e ""/4 |p( 5/4 _ .^/ 4)|2 ^ 

TT 

Even in the presence of the term in the dispersion 
relation, the power spectrum is scale-invariant in the 
case of exact de Sitter inflation. Since we have F —>■ 1 
as a; —>■ 00 , we recover the result of generalized G- 
inflation [29] in the limit e —>■ 0. For a; <C 1 we have 


N 

= No, 




(79) 

P = 

= - 

i3iV2 

hl)^M|l/ol 

< 0. 

(80) 

Since Ai > 

1 and 1 

'MY -- 

= 2/iVo2 > 0, 

we see 

that 

t/s > 0 and 

ns 

> 0. 

We also 

see that 



Fs 

0 

z 

» 


1 

_ 1 


Ml 

Ai - 


_ M|iV2 

3(2a-f 1) 
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FIG. 2: The background evolution of (a) the Hubble param¬ 
eter H and (b) the lapse function N around the Genesis-de 
Sitter transition. 


FIG. 3: (a) The sound speed squared, J-s/Qs, and (b) the 
coefficient of (divided by Qs) around the Genesis-de Sitter 
transition. 


and hence it is easy to satisfy J'g > 0 during the Genesis 
phase by choosing the parameters appropriately. 

A numerical example of the Genesis-de Sitter transi¬ 
tion is illustrated in Figs. 2 and 3. Our numerical calcu¬ 
lation was performed as follows: we solve the evolution 
equations V — 0 and df/dt = 0 with initial data {H,N) 
satisfying £ — 0, and confirm that the constraint £1 = 0 
is satisfied at each time step. In the numerical calcula¬ 
tion, the parameters are given by Mp\ = M 2 = M 3 = 1 , 
a = 1, Ai = 1 -I- 10“^, Nq = 1, and 7 = 10. The function 
f{t) is taken to be 


/ = 


A 

2 


t 


ln (2 cosh(st)) 
s 


+ /i, 


(82) 


with /o = —10“^, /i = 10, and s = 2 x 10“^. The 
background evolution is shown in Fig. 2. The evolution 
of the sound speed squared, Ts/Qsi and the coefficient 
of in the dispersion relation is shown in Fig. 3. As 
pointed out in Ref. [26], flips the sign at the transi¬ 
tion. The sound speed squared is positive except in this 
finite period. During the Genesis and subsequent de Sit¬ 
ter phases we have t /5 > 0 and > Oj and therefore we 
may conclude that this model is stable. 

Although we have thus obtained the stable example of 
the Genesis-de Sitter transition, the simple example (78) 
is not completely satisfactory if one would want successful 


gravitational reheating. Indeed, the condition (45) im¬ 
plies that X := {Ne/No)"^ < 1, but to — 4 = —2a;/(l —x) < 
0 for such X. This problem can be evaded easily by the 
following small deformation of 02 : 


1 1 + 5A^N^ 


(83) 


where A is a parameter smaller than 1/5. The condi¬ 
tion (45) now reads (1 — x)(x — 5A^) > 0, i.e., 5A^ < 
X < 1, while 


2(A -b x)(A — x) 
(1 — x)(x — 5A2) 


(84) 


is positive for 5A^ < x < A. The stability condi¬ 
tion further restricts the allowed ranges of x and A. 
The necessary condition for stability is Nf.a 2 /{N"^a 2 )' < 
0 [see Eq. (77)]. This translates to 1 -b 5A^ — 
Vl - 5A2 -b 25A4 < X < A < (4-yTT)/5 ~ 0.137, lead¬ 
ing to TO < 24/5 = 4.8. Note that the small deformation 
of 02 with A < 0.1 does not change the background and 
perturbation dynamics of the Genesis and inflationary 
phases. 

To illustrate the final stage of inflation, let us take 

i/(a+i) 

, (85) 


/ = <! + 77 


t + 


In (2 cosh(s'i)) 
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FIG. 4: The background evolution of (a) the Hubble parame¬ 
ter H and (b) the lapse function N around the end of inflation. 

where the origin of time is shifted so that the end of in¬ 
flation is given by t ~ 0. In the numerical plots pre¬ 
sented in Figs. 4 and 5, the parameters are given by 
s' = 10“^, w = 6, and A = 0.05, while the other pa¬ 
rameters are taken to be the same as the previous exam¬ 
ple of the Genesis-de Sitter transition. It is found that 
m ^ 4.5 > 4. Again, we see that < 0 in the finite pe¬ 
riod around the transition. However, Qs and T-Ls remain 
positive all through the inflation and subsequent stages. 

VI. DISCUSSION AND CONCLUSION 

In this paper, we have introduced a generic description 
of Galilean Genesis in terms of the ADM Lagrangian and 
constructed a concrete realization of inflation preceded 
by Galilean Genesis, i.e., the scenario in which the uni¬ 
verse starts from Minkowski spacetime in the asymptotic 
past and is connected smoothly to the inflationary phase 
followed by the graceful exit. Our model utilizes the re¬ 
cent extension of the Horndeski theory, which has the 
same number of propagating degrees of freedom as the 
Horndeski theory and thus can avoid Ostrogradski insta¬ 
bilities. This approach allows us to cover the background 
and perturbation evolution in all the three phases with 
the same single Lagrangian, as opposed to the effective 
field theory approach. In our scenario, the sound speed 
squared during the transition from the Genesis phase to 



FIG. 5: (a) The sound speed squared, J-s/Qs, and (b) the 
coefficient of k'^ (divided by Qs) around the end of inflation. 

inflation becomes negative for a short period. However, 
thanks to the nonlinear dispersion relation arising from 
the fourth-order derivative term in the quadratic action, 
modes with higher momenta are stable and the growth 
rate of perturbations with smaller momenta is finite and 
under control. It should also be noted that the sound 
speed of the primordial perturbations can be smaller than 
unity by choosing the parameter of the model appropri¬ 
ately. 

Although we have constructed our inflation model in 
order to resolve the initial singularity and possible trans- 
Planckian problems by incorporating Galilean Genesis 
phase before inflation, we could make use of our model 
to realize the original Galilean Genesis scenario, which is 
an alternative to inflationary cosmology, simply by taking 
vanishingly short period of inflation there. As discussed 
in the Appendix, the sound speed squared becomes nega¬ 
tive at the transition also in this case, but the instabilities 
are relevant only for small k modes thanks to the term 
in the dispersion relation. Thus, the transition from the 
Genesis phase to the reheating stage is described in a 
healthy and controllable manner. 

In fact, it would be fair to say that such a cosmol¬ 
ogy works quite well among the proposed alternatives to 
inflation, because, in contrast with the bouncing cosmol¬ 
ogy, in which all the would-be decaying modes in the 
expanding universe such as vector fluctuations and spa¬ 
tial anisotropy severely increase in an undesirable man- 
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pendix, we will go back to the original motivation of 
Galilean Genesis and study how we can match smoothly 
the Genesis phase to the reheating phase. Our approach 
based on the ADM Lagrangian is quite useful in analyz¬ 
ing such a situation as well. 


FIG. 6: The background evolution of (a) the Hubble param¬ 
eter H and (b) the lapse function N around the Genesis¬ 
reheating transition. 



ner, the Genesis solution is an attractor and generation 
of nearly scale-invariant curvature perturbation is also 
possible with an appropriate choice of model parameters 
[18]. Since no first-order tensor perturbation is generated 
in this type of scenarios, detection of tensor perturbation 
with its amplitude larger than 10“^*^ would be a smoking 
gun of inflation. 
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Appendix A: Matching Genesis to the Reheating 
Phase 

In the main text, we consider the scenario in which 
Galilean Genesis is followed by inflation. In this ap- 


FIG. 7: The sound speed squared Ts/Gs (a) and the coeffi¬ 
cient of fc'* (divided by Qs) (b) around the Genesis-reheating 
transition. 


It is now obvious that by taking 


/ 


|t|, for t < 0 

ti/(“+i), for t > 0 ’ 


(Al) 


and gluing the two functions smoothly at around t = 0, 
one can describe the Genesis-reheating transition. As a 
concrete example, we glue / ~ 0.1(—t) and / ~ 
smoothly at around t = 0 and perform a numerical calcu¬ 
lation as shown in Figs. 6 and 7. The other parameters 
are the same as those taken in the main text. As is 
expected, the numerical result here is much the same as 
the case where a duration of the intermediate inflationary 
phase is taken to be very short. In particular, becomes 
negative at the Genesis-reheating transition. The model 
is nevertheless stable since the conditions t/s > 0 and 
TLs > 0 remain satisfied. 
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